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We use the configuration interaction method and many-body perturbation theory to perform 
accurate calculations of energy levels, transition amplitudes, and lifetimes of low-lying states of 
barium and radium. Calculations for radium are needed for the planning of measurements of parity 
and time invariance violating effects which are strongly enhanced in this atom. Calculations for 
barium are used to control the accuracy of the calculations. 
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I. INTRODUCTION 

Measurements of the effects of parity (P) and time (T) 
invariance violation in atoms are an effective means to 
search for new physics beyond the standard model (see, 
e.g., Ref. pj). The radium atom is a good candidate 
for this search due to the sizeable enhancement of the P- 
and T-odd effects arising from the nuclear and electronic 
structures 0, Q . Preparations for measurements are cur- 
rently in progress at Argonne National Laboratory Q| 
andKVIpl. 

In our previous work [3j we performed a detailed study 
of those P- and T-odd effects in radium that are en- 
hanced due to the close proximity of states of opposite 
parity. Estimates of the lifetimes of all low states of ra- 
dium were also presented. 

A detailed knowledge of the energy levels and tran- 
sition amplitudes of radium is extremely important at 
the first stages of the experiment when the trapping and 
cooling scheme is developed and tested. Energy levels of 
radium presented in Moore's book [(| are based on works 
by Rasmussen Q and Russell conducted as early as 
1934. The first work presents measurements of tran- 
sition frequencies while the second work |8|] corrects the 
interpretation of these measurements. To the best of our 
knowledge no further measurements were performed for 
radium since that time. There is some concern ignited 
by recent calculations by Bieroh et al. [|| that the po- 
sitions of the energy levels of the 7s6d configuration of 
radium might be lower then those presented in Moore's 
book |(J. This, if true, can totally destroy the cooling 
scheme adopted by the Argonne group Q . 

The actual position of the 7s6d 3 D2 energy level is 
also important for the enhancement of the P- and T-odd 
effects considered in Refs. ,3, EJ- This enhancement is 
due to the very small energy interval (~ 5 cm^ 1 ) between 
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states of opposite parity 7s7p 3 P\ and 7s6d 3 £>2- Any 
significant change in the position of either state would 
also destroy the enhancement. 

We haven't performed accurate calculations of radium 
ener gy levels before. However, calculations for bar- 
ium [3, EI; which has a similar electron structure, show 
that such calculations are possible. The theoretical un- 
certainty cannot be as low as 5 cm -1 needed to confirm 
the strong P-odd enhancement due to the small energy 
interval between states of opposite parity. However, it 
can be small enough to address any concern about the 
experimental numbers. 

In the present work we perform accurate calculations 
of the energy levels and El and E2 transition ampli- 
tudes for low states of radium and barium. We use the 
V N ~ 2 approximation (see, e.g., Ref. |ll]). Relativistic 
Hartree-Fock calculations are carried out for a doubly 
ionized ion with both valence electrons removed. The 
self-consistent potential of the core (the V N ~ 2 potential) 
is used to construct the effective Hamiltonian for the 
configuration interaction (CI) treatment of the valence 
electrons. Core-valence correlations are also included by 
means of many-body perturbation theory (MBPT). Cer- 
tain classes of dominating higher-order diagrams for the 
core-valence correlation are included in all orders in the 
Coulomb interaction. 

The barium and radium atoms have similar electron 
structure, but more experimental data is available for 
barium. Parallel calculations for these atoms provide a 
control of the accuracy. The resulting accuracy for the 
energies of barium and radium is a fraction of a percent or 
better for removal energies and for the intervals between 
ground and low-lying states. There is also very good 
agreement between experimental and calculated lifetimes 
of several states of barium. This means that the accuracy 
of the results for radium should also be very high. 

In the end we see no reason for not trusting the exper- 
imental energies of radium presented in Moore's book. 
Therefore, the trapping and cooling of radium should 
work as planned. 
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II. CALCULATIONS 

We use the combined configuration interaction and 
many-body perturbation theory method (CI+MBPT, 
[l^ ) and the V N ~ 2 approximation (see ^jj) to perform 
the calculations. Like in the standard CI method, the 
Schrodinger equation for the wave function of two va- 
lence electrons is written in matrix form 

(H cS - E)V = 0. (1) 

^ is expressed as an expansion over single-determinant 
two-electron wave functions 

* = ^c l $,(r 1 ,r 2 ). (2) 

i 

$i are constructed from the single-electron valence basis 
states calculated in the V N ~ 2 potential. E in Eq. (JTJ 
is the valence removal energy (energy needed to remove 
two valence electrons from the atom). 
The effective Hamiltonian has the form 

H cS = h 1 (r 1 ) + h 1 (r 2 ) + h 2 (r 1 ,r 2 ), (3) 
hi(ri) is the one-electron part of the Hamiltonian 

Ze 2 

hi = ca ■ p + (f3 - l)mc 2 - + V N ~ 2 + E 1 . (4) 

r 

Si is the correlation potential which represents the cor- 
relation interaction of a valence electron with the core. 
It is essentially the same as for atoms with one valence 
electron (see, e.g., [Til IT3. fT^ ^>. 

hi is the two-electron part of the Hamiltonian 

e 2 

&2=| , +£ 2 (ri,r 2 ), (5) 

|ri -r 2 | 

£2 is the two-electron part of core- valence correlations. It 
represents screening of the Coulomb interaction between 
valence electrons by core electrons. 

The terms £1 and £ 2 can be calculated using MBPT. 
The expansion starts from the second order and inclu- 
sion of the second-order core- valence correlations into the 
effective CI Hamiltonian is very important for obtain- 
ing good agreement with experiment (see, e.g., 0>E])- 
However, as demonstrated in Ref. [ll]], inclusion of 
higher-order core- valence correlations leads to further sig- 
nificant improvement of the results. 

It is convenient to start the calculations from a pos- 
itive ion (Ba + or Ra + ) with one valence electron. The 
equation for a single-electron valence state v has the form 

{hi +ti- e v )ip v = 0, (6) 

where hi is given by Eq. (JIJ. Solving this equation pro- 
duces so-called Brueckner orbitals ip v and energies e v in 
which correlations with the core are included by means 
of the correlation potential Si. Comparing e v with the 



TABLE I: Energy levels of Ba + and Ra + in different approx- 
imations. Energies are given in cm -1 with respect to the 
continuum, minus sign is omitted. A c = E(exp) — E(HF), 
A = E(exp) - E(E (oo) ). 



State 


Exp. [6] 


HF 


A c 




£(00) 


A 






Barium 








6«l/2 


80687 


75339 


5348 


82318 


80816 


-129 


5d 3 / 2 


75813 


68139 


7674 


77224 


76345 


-532 


5d 5/2 


75012 


67665 


7347 


76286 


75507 


-495 


6P1/2 


60425 


57265 


3160 


61180 


60603 


-178 


6P3/2 


58734 


55873 


2861 


59388 


58879 


-145 




Radium 








7si /2 


81842 


75898 


5944 


83826 


81988 


-146 


6^3/2 


69758 


62356 


7402 


71123 


70099 


-341 


6d 5 /2 


68099 


61592 


6507 


69101 


68392 


-293 


7P1/2 


60491 


56878 


3613 


61386 


60702 


-211 


7P3/2 


55633 


52906 


2727 


56245 


55753 


-120 



experimental spectrum of the positive ion is a way to 
study different approximations for Ei. Table [I] compares 
energy levels of Ba + and Ra + , calculated in different ap- 
proximations, with experiment. Hartree-Fock (HF) en- 
ergies correspond to Ei = in Eq. ©. E' 2 ' is the 
correlation potential calculated in the second order of 
MBPT. £v°°) is the correlation potential in which two 
classes of higher-order diagrams are included in all or- 
ders. These are the screening of the Coulomb interaction 
and the hole-particle interaction. This is done in exactly 
the same way as in our previous works for atoms with one 
valence electron (see, e.g., 0,0,0,^1)- O ne can see 
from the table that inclusion of core-valence correlations 
is very important for obtaining accurate results and in- 
clusion of higher-order core-valence correlations leads to 
further significant improvement. The energies obtained 
with £(°°) are within 0.2-0.3% of the experimental values 
for s and p states, and less accurate for d states. 

The column marked A c in Tabled lists correlation en- 
ergies (the difference between Hartree-Fock and experi- 
mental energies) of valence states of Ba + and Ra + . One 
can see that the largest correlation energy is in d states. 
This means that these states are more sensitive to the 
treatment of the correlations and generally are harder to 
calculate to high accuracy. This is why the accuracy for 
d states of Ba + and Ra + is not as good as for s and p 
states. Large core-valence correlations for d states also 
manifest themselves in the energies of two-electron con- 
figurations containing d-electrons (e.g., 6s5d configura- 
tions of Ba and 7s6d configurations of Ra, see discussion 
below) . 

Calculations for positive ions give us a very good ap- 
proximation for the Ei operator in the Hamiltonian @ 
for a two-electron system. However, we also need to cal- 
culate the two-electron operator E 2 . We calculate it in 
the second-order of MBPT. Formally, the MBPT expan- 
sion for £1 and £ 2 goes over the same orders of pertur- 
bation theory. However, numerical results show that an 
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accurate treatment of Si is usually more important than 
that of £2 ■ Although inclusion of the higher-order corre- 
lations into £2 may lead to further improvement of the 
results, we leave this at the moment for future work. 



A. Energies of barium and radium 

Theoretical and experimental energies of neutral bar- 
ium and radium are presented in Table ITU Experimental 
values are taken from Moore's tables |g. We present two- 
electron removal energies (in a.u.) for the ground state 
(^Sq) of both atoms. The experimental value is the sum 
of the ionization potential of the neutral atom and its 
positive ion. Energies of excited states are given in cm -1 
with respect to the ground state. 

The CI column in Table |H] corresponds to the stan- 
dard configuration interaction method [£1,2 = in Eqs. 
Q and JSJ]. It takes into account correlations between 
valence electrons but neglects correlations between core 
and valence electrons. We use B-splines to construct the 
basis of single-electron states. 50 B-splines are calcu- 
lated in a cavity of radius 40as, where as is the Bohr 
radius. Eigenstates of the Hartree-Fock Hamiltonian are 
constructed from these B-splines and the 14 lowest states 
above the core in each of the s,p 1/2 , P3/2, 4/2 > d 5 / 2 , h/2, 
and waves are used in CI calculations. The uncer- 
tainty due to incompleteness of the basis is very low. It 
is < 10cm -1 for s 2 and sp configurations and < 50cm -1 
for sd configurations. 

The next column (A c ) lists the difference between ex- 
perimental and CI energies. This difference is mostly due 
to core- valence correlations. There are also contributions 
to A c due to the Breit interaction, radiative corrections, 
incompleteness of the basis for valence-valence correla- 
tions, etc. However, all these contributions are small. 

We include core-valence correlations by introducing 
operators £1 and £2 into the effective CI Hamiltonian 
[see Eqs. (0J and JSJ]. Comparison of the A c values for 
neutral Ba and Ra presented in Table [H] with the cor- 
relation energies (A c ) for positive ions (Table UJ reveals 
that core- valence correlations have a larger effect on the 
energies of positive ions than on neutral atoms. This is 
due to a cancellation of contributions from £1 and £2. 
The Hamiltonian for a positive ion © has only £1, while 
the Hamiltonian for a neutral atom 1151) has both. On 
the other hand, in the V N ~ 2 approximation used in the 
present work, the £1 operator for a neutral atom is the 
same as for a positive ion. 

The cancellation between the two types of core- valence 
correlations (£1 and £2) has an effect on the accuracy 
of calculations. The accuracy is poorer when the can- 
cellation is stronger. It is easy to see that the strongest 
cancellation takes place for sd configurations of barium 
and radium. Indeed, correlation corrections to the en- 
ergies of d states of Ba + and Ra + are about two times 
larger than those for s and p states (see Table UJ. How- 
ever, corrections to the energies of sd configurations of 



neutral barium and radium are about the same or even 
smaller than for s 2 and sp configurations. 

The column in Table |H] marked by £^ 2 ) lists results 
obtained with both £1 and £2 calculated in the second 
order of MBPT. We use the same B-splines to calculate 
£ as for the CI calculations. However, we use 45 out 
of 50 eigenfunctions and go up to I = 5 in the partial 
wave expansion. Inclusion of second-order £ leads to 
significant improvement of the results. The remaining 
deviation from experiment is just a small fraction of the 
total core- valence correction A c . However, we do further 
steps in trying to improve the results. We replace the 
second-order £1 with the all-order operator £^ 00 ' ) . We 
use the Feynman diagram technique as described in our 
earlier papers [3 [H [3 to calculate E^* . The results 
are presented in column £^°°' ) of Table [H] As one can 
see, inclusion of higher-order correlations into £1 leads 
to significant improvement of the removal energies but 
not of the energy intervals (see also Ref. [HI). There 
are at least two reasons for this. First, the change in £1 
operator between the positive ion and neutral atom, and 
second, higher orders in £2. £ is an energy-dependent 
operator: £ = £(e). It should be calculated at the energy 
of the state for which it is to be used. For example, £ s 
for the 6s state of Ba+ should be calculated at e = e(6s), 
etc. Using exactly the same £1 operator for the positive 
ion and neutral atom corresponds to an approximation 
in which the energy parameter for £ is chosen assuming 
that two-electron energy of a neutral atom is equal to the 
sum of the two single-electron energies of a positive ion. 
This approximation is too rough and some adjustment in 
the energy parameter is needed. An accurate adjustment 
is ambiguous. For example, the s-wave £ for the s 2 and 
sp configurations are not the same since these configu- 
rations have different energies. Moreover, £ operators 
for different states of the same configuration are not the 
same since different states have different energies. 

In the present paper we use a simpler way of adjusting 
the value of the £ operator, leaving an accurate treat- 
ment of its energy dependence for future work. We scale 
the single-electron part of the operator £1 while leav- 
ing the two-electron part £2 unchanged. Numerous tests 
show that any reasonable change in £2 does not lead to a 
significant change in the spectra of Ba or Ra. Therefore, 
we scale £1 to fit known energies of Ba, Ba + , and Ra + 
and use this scaling to calculate energies of Ra. 

The most straightforward way to scale the energy lev- 
els of Ra would be to perform an accurate fitting of the 
energy levels of Ba and use the same scaling parameters 
to do calculations for Ra. However, this method does not 
take into account the real difference in electron structure 
of the atoms. The ordering of the energy levels of Ba 
and Ra are different. States of the sd configuration lie 
below the sp configuration for Ba and above the sp con- 
figuration for Ra. Actually, there are more similarities 
between the neutral atom and its positive ion than be- 
tween neutral Ba and Ra. We can use these similarities 
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TABLE II: Ground state ( 1 So) removal energies (a.u.) and excitation energies (cm x ) of low states of barium and radium. 
A c = E(exp) - E(CI), A = E(exp) - E(/ 2 E<° o) ). 



State 


Exp. [6] 


CI 




£(2) 


£(°o) 


/i£ (oc) 




A 


Other 












Barium 












6s 2 


'Sn 


-0.55915 


-0.52358 




-0.56996 


-0.55903 


-0.55799 


-0.55915 






6s5d 


3 D 1 


9034 


11585 


-2551 


8956 


8425 


8730 


9040 


-6 






3 D 2 


9216 


11662 


-2446 


9165 


8611 


8910 


9217 


-1 






3 D 3 


9597 


11835 


-2228 


9609 


8999 


9283 


9582 


14 






l D 2 


11395 


12833 


-1438 


11733 


11020 


11323 


11627 


-232 




6s6p 


3 Po 


12266 


9947 


2319 


13088 


12377 


12400 


12270 


-4 






3 Pi 


12637 


10278 


2359 


13466 


12740 


12754 


12638 


-1 






3 P 2 


13515 


11019 


2496 


14374 


13611 


13596 


13518 


-3 






l Pi 


18060 


16919 


1141 


18631 


17778 


17832 


17834 


227 




5cZ6p 


3 F 2 


22065 


22018 


47 


22765 


21502 


21828 


22041 


23 






3 F 3 


22947 


22573 


376 


23732 


22403 


22698 


22926 


21 






3 F 4 


23757 


23172 


585 


24624 
Radium 


23230 


23500 


23746 


11 




7s 2 


'So 


-0.56690 


-0.52546 




-0.58071 


-0.56687 


-0.56567 


-0.56695 




-0.57979 a 


7s7p 


3 Po 


13078 


10380 


2698 


14202 


13277 


13293 


13132 


-53 


14268° 




3 Pi 


13999 


11240 


2759 


15118 


14161 


14166 


14027 


-28 


15159" 




3 P2 


16689 


13473 


3216 


17879 


16813 


16764 


16711 


-22 


17937° 


7s6d 


3 Dl 


13716 


15231 


-1515 


14043 


13342 


13423 


13727 


-11 


14012" 




3 D 2 


13994 


15284 


-1290 


14371 


13612 


13683 


13980 


14 


14465° 
12958 6 




3 D 3 


14707 


15461 


-754 


15254 


14323 


14364 


14642 


65 


15921° 




l D 2 


17081 


16798 


283 


18052 


17007 


17060 


17333 


-252 




7s7p 


'Pi 


20716 


18686 


2030 


21547 


20487 


20459 


20450 


266 


21663° 
20835° 


7s8s 


2 Si 


26754 


24030 


2724 


27643 


26673 


26571 


26669 


85 




6d7p 


3 F 2 


28038 


26328 


1710 


29425 


27736 


27833 


28001 


37 






3 F 3 


30118 


27713 


2405 


31745 


29848 


29891 


30077 


41 






3 F 4 


32368 


29383 


2985 


34195 


32134 


32129 


32370 


-2 





a Dzuba et al. fTjI . 
6 Bierori et al. 0|. 



to construct a fitting procedure which takes into account 
the difference between Ba and Ra. 

First, we scale Si to fit the energies of Ba + and Ra + . 
Fitting coefficients are presented in Table lifTl They are 
slightly different for d-states of Ba + and Ra + (0.94 for 
d-states of Ba + and 0.96 d-states of Ra + ). This is be- 
cause the 5d states of Ba + are closer to the core and 
the correlation correction is larger. Our £(°°) operator 
is less accurate for d states than for s and p states and 
the larger correlation correction leads to a noticeable loss 
in accuracy. However, with the value of Si reduced by 
only six or less percent, energy levels of Ba + and Ra + 
are fitted exactly. 

Then we use the same scaling of Si to calculate the 
energy levels of neutral Ba and Ra. The results are pre- 
sented in Table ITTl under the /iE^°°' > mark. There are 
two important things to note. The first is the significant 
improvement in the agreement with experiment. The 
second is the remarkable similarity between Ba and Ra 
which was never that good for any other approximation 
used so far. Now all 3 D states of both atoms are about 
300 cm -1 below the experimental values while all 3 P 



states are about 100-200 cm -1 above the experimental 
values. This is enough to indicate that the experimen- 
tal energies of Ra are correct or that at least there is no 
reason to believe otherwise. 

However, we do one more step. We change the scal- 
ing parameters by fitting to the energy levels of neutral 
barium. The new values are presented in Table ITTT1 This 
change in the scaling parameters accounts for the energy 
dependence of the Si operator discussed above. Then 
the scaling parameters for Ra are calculated using the 
formula 



/i(Ra) = MBa) 



/i(Ra+) 
/i(Ba+)' 



In other words, barium scaling parameters are cor- 
rected using the difference in fitting for Ba + and Ra + . 
The new fitting parameters for Ra are also presented in 
Table imi 

The results of calculations for Ba and Ra with the new 
fitting parameters are presented in Table [H] under the 
f2^i°^ mark. One can see that the 1 S, 3 D, and 3 P states 
of Ba are fitted almost exactly. The 1 D and 1 P states are 
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TABLE III: Fitting factors / for rescaling of E (oo) to repro- 
duce experimental energies of Ba + and Ra + 



Atom 


«l/2 


Pl/2 


P>3/2 


^3/2 


^5/2 


Ba II 


0.9777 


0.95 


0.96 


0.94 


0.94 


Ba I 


1.0032 


1.046 


1.046 


0.9164 


0.9164 


Ra II 


0.9777 


0.95 


0.96 


0.96 


0.96 


Ra I 


1.0032 


1.046 


1.046 


0.9359 


0.9359 



less accurate because strictly speaking one cannot use the 
same E for 3 P and 1 P states and for 3 D and 1 D states 
due to the difference in energies of these states. 

Calculations for Ra with the new scaling parameters 
reduce the deviation of the theoretical values from ex- 
periment to about 50 cm -1 or less for all low states. Let 
us stress once more that no knowledge of the Ra spec- 
trum was used to do the fitting. Values of the scaling 
parameters were found by fitting the spectra of Ba, Ba + , 
and Ra + . The good agreement of the final numbers with 
experiment leaves no room for any claim that the experi- 
mental values might be incorrect. As discussed by Russell 
8] , the difference between two possible ways of interpre- 
tation of the experimental data is 627.66 cm' 1 . This is 
much larger than the difference between our calculated 
energies and the experimental energies of Ra. 

In the last column of Table [H] we present the results 
of our previous calculations of energy levels of Ra 0J 
together with the results of recent calculations for Ra by 
Bierofi et al. [9j. Our previous calculations were very 
similar to those presented in the table in the E^ 2 ) col- 
umn. They were also obtained in the V N ~ 2 approxi- 
mations with the second-order E. However, the basis of 
single-electron states was different. The difference be- 
tween the present E^ 2 ) results and the results of Ref. 0] 
can serve as an upper limit on the uncertainty due to in- 
completeness of the single-electron basis for the valence 
states. The real uncertainty of the present calculations 
is several times smaller due to better saturation of the 
basis. 

Calculations of Ref. were performed by means of 
the multi-configuration Dirac-Hartree-Fock method. The 
authors use the results to claim that experimental en- 
ergies of the 7s6d configuration of Ra might be incor- 
rect. Indeed, their calculated value for the 3 i?2 state is 
1034 cm -1 below the experimental value. On the other 
hand, the deviation from experiment of the only other 
calculated energy level, the energy of the 1 P\ state, is 
only 119 cm -1 . No other energy levels of Ra were calcu- 
lated and no calculations for other two-electron systems 
were used to control the accuracy. Therefore, it is hard 
to make any judgement about the quality of these calcu- 
lations. However, let us remind the reader that calcula- 
tions for sd configurations are more difficult than for sp 
configurations due to the larger correlation interaction 
of the d electron with the core and the stronger cancel- 
lation between Ei and E2 terms (see discussion above). 
Therefore, the accuracy obtained for the 1 P\ state cannot 



serve as a guide for the accuracy for the state. Apart 
from that, good agreement with experiment for just one 
number cannot rule out a fortunate coincidence. 



B. Transition amplitudes 

The leading contribution to the amplitude of a transi- 
tion between states v and w of Ba or Ra is given by 

A vw = (* w \f\* v ), (7) 

where ^ w and ^> v are the solutions of Eq. Q and / 
is the operator of the external field. This expression 
doesn't take into account the effect of the external field 
on the atomic core. This effect, which is known as core 
polarization, is very important and can change the am- 
plitude significantly. It can be included by means of the 
time-dependent Hartree-Fock method (TDHF) which is 
equivalent to the well known random-phase approxima- 
tion (RPA) method. 

Every single-electron core function is presented in the 
RPA approximation as i/> a +(5V>a, where ip a is the Hartree- 
Fock wave function of the core state a calculated in the 
V N ~ 2 potential; 6ip a is the correction due to the external 
field. The corrections to all core states are found self- 
consistently by solving Hartree-Fock-like equations 

{Hq - e^Slpa = -flp a - 5V core ll>a,, (8) 

where Hq is the Hartree-Fock Hamiltonian, / is the oper- 
ator of the external field, and 8V core is the correction to 
the self-consistent potential of the core due to the effect 
of the external field. Note that in our case V core = V N ~ 2 . 
The SV core term is calculated using the Sip corrections to 
all core states. The final expression for the transition 
amplitude has the form 

A vw = (V w \f + 6V core \y v ). (9) 

Amplitudes of electric dipole transitions (El) between 
low states of barium and radium calculated in different 
approximations are presented in Tabic IIVI Core polar- 
ization is included everywhere since it is known to be 
an important effect. We study only the effect of core- 
valence correlations on the amplitude. As with energies, 
inclusion of core-valence correlations have a significant 
effect on the amplitudes. On the other hand, amplitudes 
calculated with E^ 2 ' and E^ 00 ) are not very different. 

We also present in Table IIVI the results of our previ- 
ous calculations for El transition amplitudes Q. In spite 
of the very simple approximation for the wave functions 
used in the previous work, the agreement for the ampli- 
tudes is generally remarkably good. The exception is the 
amplitudes which involve a change of spin. These ampli- 
tudes are larger in the present calculations than in our 
previous work. The reason is the underestimation of rela- 
tivistic effects for the 5d state of Ba and 6d state of Ra in 
Ref. [3 • The electric dipole transitions between states of 
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TABLE IV: El-transition amplitudes for Ba and Ra in differ- 
ent approximations (|(i||dz||j)|ao)- 



Transition 


CI 


t (2) 




Other [31 








Barium 






3 Po 


3 D! 


2.6185 


2.3149 


2.3045 


2.3121 


3 Pi 


'So 


n 3903 


0.5281 


0.5240 


0.4537 


3 P! 


3 Dr 


9 9899 


2.0104 


9 0096 

_■ . \J\J 


2.0108 


3 Pl 


3 D 2 


3 8806 


3.4309 


3.4128 


3.4425 


3 Pl 


'D 2 


0.2979 


0.4675 


0.4999 


0.1610 


3 P 2 


3 Di 


0.5997 


5269 


0.5247 


0.5275 


3 P 2 


3 D 2 


9 9838 


2.0012 


1.9933 


2.024 


3 P 2 


3 D 3 


5.4285 


4.8181 


4.7805 


4.777 


3 P 2 


'D 2 


0.3321 


0.3551 


0.3402 


0.1573 


'Pi 


'So 


5.7133 


5.4235 


5.4695 


5.236 


'Pi 


3 Di 


0880 


0850 


0.0735 


0.1047 


'Pi 


3 D 2 


5935 


0.4143 


3999 


0.4827 


-1 1 




0.9919 


1.3062 
Radium 


1.1394 


1.047 


3 Po 


3 Di 


3.2996 


2.9325 


2.9521 


3.0449 


3 Pl 


'So 


0.8241 


1.2317 


1.2205 


1.0337 


3 Pl 


3 D! 


2.8836 


2.5155 


2.5366 


2.6389 


3 Pl 


3 D 2 


4.8393 


4.2931 


4.3158 


4.4399 


3 Pl 


'D 2 


0.7095 


0.7397 


0.8068 


0.0467 


3 P 2 


3 D! 


0.7799 


0.6714 


0.6781 


0.7166 


3 P 2 


3 D 2 


2.9438 


2.5357 


2.5615 


2.7283 


3 P 2 


3 D 3 


6.9465 


6.2626 


6.2541 


6.3728 


3 P 2 


'D 2 


0.4285 


0.5885 


0.5344 


0.1499 


'Pi 


'So 


5.7703 


5.3652 


5.4821 


5.4797 


'Pi 


3 D! 


0.3736 


0.4381 


0.4255 


0.4441 


'Pi 


3 D 2 


0.6162 


0.3965 


0.3591 


1.188 


^1 


'D 2 


2.9923 


3.3103 


3.1379 


2.4053 



different spin are forbidden in the non-relativistic limit. 
Therefore, larger amplitudes means larger relativistic ef- 
fects. Since we don't have experimental values for the 
amplitudes, we can use fine structure intervals instead 
to see how well relativistic effects are treated in different 
calculations. One can see from the data given in Table I 
of Ref. Q that the fine structure intervals between the 
3 -Di,2,3 states of Ba and Ra are about two times smaller 
than the experimental values. In contrast, all fine struc- 
ture intervals of the present calculations are very close 
to experiment (see Table [HJ). Therefore, we expect the 
corresponding amplitudes to be more accurate. 

C. Lifetimes of low- lying states of barium and 
radium 

The lifetime of the atomic state i expressed in seconds 
is given by 

n =2.4189 x 10- 17 /Y^T tJ , (10) 

3 



where Ty is the probability of a transition from state i to 
a lower state j (in atomic units), the numerical factor is 
to convert atomic units to seconds, and summation goes 
over all states j that have energies lower than the energy 
of state i. 

In the present paper we consider only electric dipolc 
(El) and electric quadrupole (E2) transitions. The prob- 
ability of the El transition from state i to a lower state 
j is (atomic units) 



where Uij = 6j — tj , Ajj is the amplitude of the transition 
(reduced matrix element of the electric dipole operator), 
and Ji is the value of the total angular momentum of the 
state i. The probability of the E2 transition is (atomic 
units) 



Lifetimes of low- lying states of barium and radium cal- 
culated using transition amplitudes from Table IIVI and 
experimental energies are presented in Tables IVl and IVll 
The new data show systematic improvement of the agree- 
ment between theory and experiment compared to our 
previous work However, the change is small. The 
only significant change is for the 3 Pi states of barium 
and radium and 3 Z?2 state of radium. This is due to the 
change in the El amplitude of the 3 Pi — 1 So transition 
and the E2 amplitude of the 3 Z?2 — 1 5'o transition. The 
new values are more accurate due to the better treatment 
of relativistic effects (see discussion above). 
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TABLE V: Lifetimes of low-lying states of barium. Numbers in square brackets denote powers of 10. 



State 


State 


Transition 


Lower states to decay to 

u> Amplitude 
(a.u.) (a.u.) 


Probability (s 1 ) 
This work Ref. f!8l 


This work 


Lifetime 
Dzuba et al. [3] 


Other 


3 Po 
u 


3 Di 


El 


0.01454 


2.305 


3.500[ 5] 




2.86 its 


2.83 us 




3 P1 


% 


El 


0.05758 


0.5241 


3.743[ 5] 




1.23 us 


1.37 us 


1.2 us" 




3 Di 


El 


0.01642 


2.003 


1.267[ 5] 












3 D 2 


El 


0.01559 


3.413 


3.150[ 5] 












J D 2 


El 


0.00567 


0.4999 


3.257[ 2] 










3 P 2 


3 Dr 


El 


0.02041 


0.5247 


1.003[ 4] 




1.44 us 


1.41 us 






3 D 2 


El 


0.01958 


1.993 


1.278[ 5] 












3 D 3 


El 


0.01785 


4.781 


5.570[ 5] 












X D 2 


El 


0.00967 


0.3403 


4.484[ 2] 










X P1 


X S 


El 


0.08238 


5.470 


1.194[ 8] 


1.19(1)[8] 


8.35 ns 


9.1 ns 


8.37(8) ns 6 




3 Dr 


El 


0.04122 


0.0735 


2.702[ 3] 


3.1(1)[3] 










3 D 2 


El 


0.04039 


0.3993 


7.498 [ 4] 


1.1(2)[5] 










J D 2 


El 


0.03047 


1.139 


2.623[ 5] 


2.5(2) [5] 








3 D 2 


X S 


E2 


0.04199 


3.125 


1.454[-2] 




69 s 







"Reference |19|. 
'Reference |2(| . 



TABLE VI: Lifetimes of low-lying states of radium. Numbers in square brackets denote powers of 10. 



State 


State 


Lower states to decay to 
Transition uj Amplitude 
(a.u.) (a.u.) 


Probability 


This work 


Lifetime 
Dzuba et al. [3J 


Other 


3 D! 


3 Po 


El 




0.00291 


2.952 


1.529[ 3] 


654 us 


617 us 




'Da 


3 Pi 


El 




0.01404 


0.8068 


7.722 [ 3] 


129 ms 


38 ms 






3 P 2 


El 




0.00179 


0.5345 


7.973 [ 0] 








3 Pr 


% 


El 




0.06378 


1.221 


2.760[ 6] 


362 ns 


505 ns 


420 ns°,250 ns 6 




3 Di 


El 




0.00129 


2.537 


9.850[ 1] 










3 D 2 


El 




0.00002 


4.316 


1.572[-3] 








3 P 2 


3 Dx 


El 




0.01355 


0.6782 


4.897[ 3] 


5.55 us 


5.2 ps 






3 D 2 


El 




0.01228 


2.562 


5.204[ 4] 










3 D 3 


El 




0.00903 


6.254 


1.234[ 5] 








'Pi 


% 


El 




0.09439 


5.482 


1.805[ 8] 


5.53 ns 


5.5 ns 






3 Di 


El 




0.03189 


0.4256 


4.195[ 4] 










3 D 2 


El 




0.03063 


0.3592 


2.646[ 4] 










X D 2 


El 




0.01656 


3.138 


3.194[ 5] 








3 D 2 


'So 


E2 




0.06376 


5.022 


3.032[-l] 


3.3 s 


15 s 


4 s c 



"Reference ^1 . 
'Reference 25j. 
c Reference 31 ■ 
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